Interference effects between two initially independent Bose-condensed gases 
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When two initially independent Bose-condensed gases are allowed to overlap, we investigate the 
density expectation value of the whole system by using the second quantization method. In the 
presence of interatomic interaction, based on the exact expression of the density expectation value, 
it is found that there is a nonzero interference term in the density expectation value of the whole 
system. The evolution of the density expectation value is shown for different coupUng constants. 
The present work shows clearly that there is an interaction- induced coherence process between two 
initially independent condensates. 
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A. Introduction 



After the experimental realization of Bose-Einstein condensates about ten years ago , the intensive exper- 

imental and theoretical studies have shown clearly that the Bose condensate has stable coherence property 4, 5,l6(. 
After two separated condensates are allowed to overlap, the coherence property of the whole system is a very interest- 
ing problem. Interference patterns between Bose condensates can give us important information about the first-order 
coherence properties of the whole system. Thus, in the celebrated experiment by Andrews et al 7], the interference 
patterns between two separated condensates are investigated. In this experiment, the interference patterns between 
two coherently separated Bose condensates were observed clearly. In the same experiment, there is also a quite striking 
experimental result that there are high-contrast fringes even for two completely independent condensates at an initial 
time. 

In our recent theoretical paper , based on the exact expression of the density expectation value calculated from the 
many-body wavefunction of the whole system, it is shown clearly that there is an interaction-induced coherence process 
when two initially independent condensates are allowed to overlap. Due to the interaction-induced coherence process 
between two initially independent condensates, there would be interference patterns in the density expectation value 
when there is an overlapping between two initially independent condensates. In the present work, we will investigate 
the density expectation value based on the second quantization method. It is shown that the result calculated from 
the second quantization method is the same as the result based on the many-body wavefunction. The calculations 
based on the second quantization method will give us further physical picture for the interaction-induced coherence 
process. In particular, to show clearly the interaction-induced coherence process, we calculate the density expectation 
value for different coupling constants in this work. It is shown clearly that increasing the interatomic interaction will 
enhance the coherence effect of the whole system. 

B. The exact expression of the density expectation value 

Two separated condensates can be created by a double- well trap (See for example the experiments of Refs 0,0). 
When the tunneling effect between two wells can be omitted, one can obtain two completely independent condensates 
with laser cooling and the following evaporative cooling in the presence of the double- well trap. Obviously, when the 
tunneling effect between two wells is obvious, the two condensates confined in the double-well trap are coherently 
separated. The interference patterns of the whole system can be investigated by switching off the double- well trap. As 
shown in Fig. 1(a), there is no overlapping between two initially independent condensates. After the confined potential 
is switched off, two condensates will expand freely and result in the overlapping. 

For two initially independent condensates comprising particle number A^i and A^2, the corresponding state is (see 
also [iSlill): 
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FIG. 1: Figure 1(a) shows two initially independent condensates, while figure 1(b) shows the overlapping of two condensates 
due to free expansion. 



where Cn is a normalization constant to assure {Ni, N2\Ni, N2) = 1. a[ (a^) ^ creation operator which creates a 
particle described by the single-particle state {(j)2) in the left (right) condensate. 
One should note that for two initially coherently-separated condensates, the state is 



\N) = 



(6t)^|0) 



(2) 



where is a creation operator which creates a particle with the single-particle state {\/Ni<pi + \^T^(f>2) /V^i + -^2- 
For two coherently-separated condensates, the density expectation value is 



ne(r,t) = {N,t\^^{r,t)^{r,t)\N,t) 



= N 



ae \^i (r, t)r + 25e x Re {^1 (r, t) (f>2 (r, t)) + Ce |(/)2 (r, t)[ 



(3) 



where Oc = Ni/N, be = \/N\N2/N and Cc = N2/N with N = Ni + N2- The second term in the above equation 
accounts for the interference effect when there is an overlapping between two condensates upon expansion. For two 
coherently-separated condensates, the interference patterns were investigated in several interesting theoretical works 
12, 13]. 

For two initially independent condensates, however, for the state given by Eq. it seems that there is no 

interference term in the density expectation value with the following simple calculation (See also 



nd{r,t) = {Ni,N2,t\^^ {r,t)^ {r,t)\Ni,N2,t) 

= {Ni,N2,t\ (a\a,\<Pi{r,tf +ala2\4>2{r,tf^ \NuN2,t) 

+2 X Re [{Ni,N2,t\ a\a2 \Ni,N2, t) 4>l (r, t) 4>2 (r, t)^ 
= N^\(t>i{v,tf +N2\4>2(jC,t)\\ 



(4) 

(5) 
(6) 
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In the above equation, the field operator ^ (r, t) is expanded as ^ (r, t) = ai^i + a202 + • • • with ai and a2 being 
the annihilation operators. One should note that, to get from there is an implicit assumption that ai and a\ 
are commutative. This holds when J (j)i (r, t) (r, ^) = 0. When ai^a^ = 0, it is easy to understand that the 

interference term (the last term in Eq. is zero in rid (i*?^)- 

In the present work, however, we will show that and 02 will become non-orthogonal in the presence of interatomic 
interaction. Thus, we will investigate the general case for J 0i05(iF = (. The operators ai and a2 can be written as 



ai 



and 



a2 



I 



^4>ldV. 



(7) 



(8) 



Here ^ (x, t) is the field operator. By using the commutation relations of the field operators (x, t) , ^ (y, t)] = 
and (x, t) , ^"f" (y, t)] = 5 (x — y), it is easy to get the commutation relation 

[ai,4]=C. (9) 

We see that ai and are not commutative any more for J ^i^JdV being a nonzero value. In this case, it is obvious 
that one can not get the result from (Q. This means that one should be very careful to get the correct density 
expectation value for two initially independent condensates. 

It is well-known that the field operator should be expanded in terms of a complete and orthogonal basis set. 
Generally speaking, the field operator ^ (r, t) can be expanded as: 

$ (r, t) = aic^i (r, t) + kc^'^ (r, t) + • • • , (10) 
where and (j)'2 are orthogonal normalization wavefunctions. Assuming that (j)'2 = ((/>2 + <^(/>i ) , based on the 



conditions / 010^2^^^ = and / 1^ dV = 1, we have = (l - \C?) 
we have 

a2 k/P* + (ai. 
It is easy to get the following commutation relations: 

[k,k] = [^^^^] = 0, [^,^^] = 

[ai,ai] = [a\,a\]=0,[ai,a\] 



and a = — Based on 



(11) 

(12) 



[/c, ai] — [k,a\] — 0. 



(13) 



Because k and a[ are commutative, it is convenient to calculate the density expectation value Ud (r, t) by using the 
operators k and a[. The exact expression of the density expectation value is 



where the coefficients are 



ad 



ad (r, t) r + 2 X Re {Ml (r, t) ^2 (r, t))+jd\^2 (r, t)\' 



^N2\iN^+i-l)\N^[l-\CfY" 'iCI 
^ iUliNi - l)\{N2-iy. 



2i 



/3d = 



Id 



N2-1 N2\{Ni+i 



)!(l-|Cl') 



i=0 



i! (i + 1)! (A/"! - 1)! {N2-i- 1)! 



^iiV2!(iVi+i)!(l-|C|') 



ICI 



2i 



j=0 



(iV2 - i - 1)! 



(14) 

(15) 
(16) 
(17) 
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In addition, the normalization constant is determined by 

N2 

^ i^N2\{N^+i)\(i-\c{t)r: 



ililNil {N2 - i)\ ^^^^ 

The above density expectation value is obtained based on the second quantization method. Although the expressions 
of the coefficients given by Eqs. (|T5j) . (|T6j) . and (|T7|) are quite different from the results calculated from the many-body 
wavefunction (See Ref. lH), we have proven that this density expectation value given by Eq. (|T^ is equal to the 
result calculated from the many-body wavefunction ^at^ats (ri, • • • , tati+ats, which satisfies the exchange symmetry 
of identical bosons in Ref. 8]. 

For two independent ideal condensates, before the overlapping of the two condensates, ( {t = 0) = 0. Based on the 
Schrodinger equation, it is easy to verify that after the double-well potential separating two condensates is removed, 
one has ^ (^) = at any further time. Thus /3d = 0^ and the density expectation value given by Eq. ([T^ is equal to 
the result given by Eq. Q. In this case, the interference term is zero in the density expectation value. 

In the presence of interatomic interaction, we now turn to investigate the evolution equations of 0i and ^2- The 
overall energy of the whole system is 

E = J dV{Ni,N2,t\ (^|!-V$t. v$ + Fext$^$ + |$^$^$$) \NuN2,t), (19) 

where Vext is the external potential and g is the coupling constant. By using the ordinary action principle and the 
above interaction energy, one can get the following coupled evolution equations for ^1 and 02- 

d^i ^ _]_5E_ 

dt Ni sri ' 

•^^^2 1 SE 

= 7V2^' ^''^ 

where 5E/5(j)\ and 5E/5(j)2 are functional derivatives. With the above coupled equations, one can understand that 
C (t) becomes nonzero after the overlapping between two condensates for g being nonzero. The reason that C (t) 
becomes nonzero is especially due to the nonlinear interaction of the whole system. 

Although generally speaking, \Cif)\ is much smaller than 1 because is an oscillation function about the space 
coordinate, nevertheless, a nonzero value of (t) will give significant contribution to the density expectation value for 
large A^i and N2. It is found that for |C| being the order of N^^ and ^ the interference term in Eq. (|T^ will play 
very important role in the density expectation value |^ . The coherence effect of two initially independent condensates 
can be shown through the value of \l3d/ocd\' For Pd/o^d = 0, there is no coherence between two condensates even there 
is an overlapping. For \l3d/oid\ = 1, the two condensates can be regarded to be fully coherent. 

We see that interatomic interaction plays an essential role in the emergence of the interference effect for two initially 
independent condensates. In the presence of interatomic interaction, the two initially independent condensates will 
become coherent after the overlapping between two condensates. This interference effect is a natural result of the 



commutation relation 



ai, a2 



which is nonzero for being nonzero in the presence of interatomic interaction. 

Generally speaking, increasing the particle number will enhance the interference effect in the density expectation 
value. Based on Eq. (|2n|), increasing the coupling constant g will have the effect of increasing ( (t). Thus, increasing 
the interaction between particles will enhance the coherence effect. To show more clearly the interaction-induced 
coherence process between two initially independent condensates, here we will investigate the density expectation 
value for different coupling constants. In the last few years, the rapid experimental advances of Feshbach resonance 
[1^ where the scattering length can be tuned from positive to negative make this sort of experiment be possible. 



C. Evolution of the density expectation value for different coupling constants 



To give a clear presentation, we consider the evolution of the density expectation value for one-dimensional case. 
At t = 0, to give a general comparison for different coupling constants, the initial wavefunctions for two independent 
condensates are assumed as 



(j)! {xi,t = 0) 



1 



: exp 



{xi - Xif 
2Al 



(21) 



FIG. 2: After two initially independent condensates are allowed to expand freely, shown is the evolution of the density 
expectation value nd{xi,T) (in unit of A^i + A^2) for different coupling constants gi. Shown in the inset of each figure is 
the relation between f3d/ad and dimensionless time r. For two ideal condensates shown in figure 2(a), we see that there is no 
interference pattern even there is an overlapping between two condensates. For the case of = 1 shown in figure 2(b), we see 
that low-contrast interference patterns begin to emerge due to the interaction- induced coherence process. In figures 2(c) and 
2(d), we see that there are high-contrast interference patterns. In particular, in figure 2(d) for gi = 20, two initially independent 
condensates can be regarded to be fully coherent near r = 1. 



h{xi,t = 0) = =exp 

In the above wavefunction, we have introduced a dimensionless variable xi = x/l with / being a length. In the present 
work, we assume that Ai = A2 = 0.5, and X2 — xi = 4.5. For these parameters, at t = 0, the two condensates are 
well separated. In the numerical calculations of the coupled equations given by Eq. (|2n|), it is useful to introduce the 
dimensionless variable r = Eit/H with Ei = h'^/2mP and dimensionless coupling constant gi = Nig/Eil. In addition, 
the particle number is assumed as A/"! = A/'2 = 1.0 x 10^. In real experiments, interatomic interaction will play very 
important role in the ground-state wavefunction of the condensates confined in the double- well trap. However, in 
principle, one can prepare the state given by Eqs. (|2TJ and by adjusting the double- well trap for different coupling 
constants. In the present theoretical work, the identical initial wavefunction for different coupling constants will be 
helpful in the comparison of the density expectation value for different coupling constants. 

With these parameters, one can get the evolution of (pi and (j)2 based on the numerical calculations of Eq. (|20|) after 
the double- well potential is switched off. From 0i and ^2, we can get and thus the density expectation value based 
on Eq. (|T^ . Shown in figure 2 is the evolution of the density expectation value for different coupling constants. It is 
shown clearly that increasing the coupling constant will enhance the coherence effect. 

D. Summary and discussion 

In summary, by using the second quantization method, we give the exact expression of the density expectation value 
when two initially independent condensates are allowed to overlap. The exact expression of the density expectation 



{xi - X2) 
2Ai 



(22) 
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value is the same as the result calculated from the many-body wavefunction 8]. In the presence of interatomic 
interaction, the wavefunctions of two condensates will become non-orthogonal after the overlapping. For large particle 
number, the non-orthogonal property of the wavefunctions plays very important role in the interference effect between 
two condensates. Through the numerical calculations of the density expectation value of the whole system for different 
coupling constants, we show the interaction-induced coherence process between two initially independent condensates. 
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